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Patterns in the weak shear 3-point correlation function 
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Abstract. We explore the scale and angular dependence of the cosmic shear three-point correlation function. 
The shear field is found to have a much more complex three-point function than the convergence field, but it 
also exhibits specific motifs which draw signatures of gravitational lensing. Exact shear patterns are inferred 
analytically for some geometrical shear triplets configurations, when simple interpretations can be derived. A 
more complete description of their geometrical properties has then been carried out from ray-tracing numerical 
simulations and simple models. These patterns reveal typical features that can be utilized for non-gaussian signal 
detection in cosmic shear surveys. We test the robustness of these properties against complex noise statistics 
and non-trivial survey topologies. From these conclusive checks, we predict that the VIRMOS-DESCART survey 
should allow a detection of a non-gaussian signal with a comfortable significance for a low matter density Universe. 
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1. Introduction 

The detections of cosmic shear signal (van Waerbeke et al. 
2000, Bacon et al. 2000, Wittman et al. 2000, Kaiser et al. 
2000, Maoli et al 2001) opened the analysis of large-scale 
mass distribution in the Universe to weak lensing surveys. 
As the surveys sizes progressively increase, the noise and 
systematics are reduced to reasonable levels, which per- 
mits to constrain cosmological parameters from the am- 
plitude and the shape of the shear two-point correlation 
function (Maoli et al. 2001, van Waerbeke et al. 2001b). 
However, although it is weakly sensitive to the cosmolog- 
ical constant fl\ (Bernardeau, van Waerbeke & Mellier 
1997), the two-point function is a degenerate combination 
of the matter density f2 m of the Universe and the ampli- 
tude of the power spectrum erg (Villumsen 1996, Jain & 
Seljak 1997). 

This degeneracy can be broken from the angular depen- 
dence of the cosmic shear amplitude (Jain & Seljak 1997), 
but it relies on the prior knowledge of the shape of the 
mass power spectrum. An alternative is to directly probe 
weak shear maps, which contain non-gaussian features 
that can be used to derive f2 m with the only assump- 
tion that initial conditions were gaussian (Bernardeau, 
van Waerbeke, Mellier 1997, van Waerbeke, Bernardeau, 
Mellier 1999) - this later assumption being eventually 
testable in the data set itself . 



So far, all the theoretical predictions regarding non- 
gaussian features are based on a critical reconstruction 
process of either a convergence map (filtered for instance 
with a top-hat window function) or an aperture mass 
map (Kaiser et al. 1994, Schneider 1996, Schneider et al. 
1998, Bernardeau & Valageas 2000). Unfortunately, the 
panoramic reconstruction of mass maps from real data 
turned out to be considerably more difficult than ex- 
pected. The masking process discussed by van Waerbeke 
et al. 2000 produces patchy surveys with a non-trivial 
topology and inhomogeneous noise. The resulting mass 
maps have poorly understood statistical properties which 
is practically difficult to handle with confidence. We there- 
fore explore another option which uses direct signatures of 
gaussian effects in shear map patterns. 

Shear pattern study is a priory more difficult because 
the third order moment of the local shear vanishes for ob- 
vious symmetry reasons. It is therefore necessary to seek 
for peculiar configurations (or geometries) of shear triplets 
for which a significant signal is expected. This is the aim 
of this paper. In the next section, we present the exact 
analytical results obtained for specific geometries, their 
interpretation and the results obtained from ray-tracing 
simulations which exhibit the complete three-point func- 
tion patterns. The detection of these features in these sim- 
ulations is presented in section 3. These results suggest a 
detection strategy in cosmic shear surveys which is put 
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forward and tested with mock catalogs in section 4. We 
checked carefully that the masks do not compromise the 
possibility of detecting the effect. Finally, we estimate the 
expected error level of such measurements in the available 
data sets. 

2. Theoretical insights 

2.1. The local shear statistical properties 

In the following we assume the shear signal is entirely 
due to cosmic shear effects, and there is no so-called B- 
mode produced by lens-lens coupling, clustering or any 
other contamination. Later in the paper (in Section 4) we 
study the impact of residual systematics. This hypothesis 
naturally leads to derive the shear field 7 from a potential 
field which can be written, in Fourier space, 



7(«0 = /d»i* «p(H.*)«(ft) 



with 



u Wi 



cos(26»0 
sin(2ft) 



(1) 



(2) 



where Hi is the Fourier transform of the local convergence 
map, the components of I being lcos(9i) and ls'm(6i). We 
also assume the validity of the small angle approxima- 
tion, making the spherical harmonic decomposition un- 
necessary (Stebbins 1996). 

The statistical properties of the variable ki can be in- 
ferred from those of the projected density contrast. These 
derivations have been studied in details in the many pa- 
pers mentioned in the introduction. Using the small an- 
gle approximation (Limber 1954, applied here in Fourier 
space, Kaiser 1992, 1998, Bernardeau 1995), they give the 
expression of the convergence power spectrum P K (l), de- 
fined as, 



(KlKv) = <5Dirac(J + P K (l) 



(3) 



as an integral over the line-of-sight of the 3D matter den- 
sity contrast power spectrum P$(k), 
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where w{x) is the lensing efficiency function (x is the ra- 
dial distance, V is the comoving angular distance). It de- 
pends on the cosmological parameters and on the source 
redshift distribution function, n(z s ): 



w (x) = 



dz s n(z s ) 



V(Xs - X) v{x) 

aV( Xs ) 



(5) 



where \ s is the radial distance to the source plane and 
when the distances are expressed in units of c/Hq. The 
phcnomcnological properties of the convergence field, or 
similarly the shear field, will be determined by both the 
shape of the efficiency function and the properties of the 
cosmic density field. 



Likewise, the convergence three-point function can be 
derived from the cosmic matter three-point function. The 
small angle approximation can indeed be used again at 
this level (Bernardeau 1995, Bernardeau et al. 1997), lead- 
ing to the expression of the Z-space three point function 

(k^k^k^) =5ui Ia c(h + h + h)B K (li,l2,h) (6) 
with 



B K (h,l 2 ,h) = 

Jo 



(X)S 5 I-,-,- ) (7) 



where B$(ki,k 2 ,k 3 ) is the 3D matter density contrast 
bispectrum: 

(<5fc 1 <5fc 2 <5fc 3 ) = fcirac(fel + k 2 + k 3 ) B s {k 1 ,k 2 , k 3 ). (8) 

The bispectrum B$ is usually parametrized as follows 

B s (k 1} k 2 ,k 3 ) = Q{k 1 ,k 2 )Ps(h)Ps(k 2 ) + sym. , (9) 

where Q is an homogeneous function of the wave vectors 
k\ and k 2 . For instance, in the quasi-linear regime it is 
given by, 

10 t fei.fc 2 , kt.k 2 i 4(fc!.fc 2 ) 2 

y 



Q(ki,k 2 ) 



(10) 



whereas in the nonlinear regime Q is often assumed to be 
a pure number (that depends however on the cosmological 
models). Several phenomenological models have been pro- 
posed to describe the behavior of the three-point matter 
correlation function, from the quasi-linear regime to the 
nonlinear regime (the so-called EPT, Colombi et al. 1997 
or its extensions with the HEPT, Scoccimarro & Frieman 
1999). In particular Scoccimarro & Couchman (2001) have 
shown that the matter three-point function could be de- 
scribed by, 



Q{ki,k 2 ) = a(ki,k 2 ) + b(ki,k 2 ) 



k\.k 2 fei.fe- 



k 2 



-c(ki,k 2 



(fc!.fc 2 ) 2 

""2 



(11) 



where a, b and c depend on scale in a known way. 

The properties of the 3D matter three-point function 
actually extend to the three-point correlation function of 
the convergence field; e.g. one expects to have 



B K (h,l 2 ,l 3 ) = Q K (h,l 2 )P K (h)P K (l 2 ) + sym. 



(12) 



To a large extent the functional dependence of the three- 
point function coefficient Q K is left unchanged (Toth, 
Hollosi & Szalay 1989, Bernardeau 1995) although its am- 
plitude is affected by projection effects. This has been in- 
vestigated in detail at the level of the convergence reduced 
skewness, which is proportional to the angular averages of 
Q K (ki, k 2 j^ Interestingly the reduced skewness, defined 
as 



«3 



(* 3 ) 



(13) 



1 In the non-linear regime one simply has S3 ~ 3 Q K to a 
very good approximation 
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depends on the efficiency function amplitude, and there- 
fore on the cosmological parameter, f2 m , but is less sen- 
sitive to the cosmological constant value or the vac- 
uum equation of state (Benabed & Bernardeau 2001). 
For sources at redshift about unity, calculations give 
(Bernardeau et al. 1997, Hui 1999, van Waerbeke et al. 
2001a), 



S3 



(14) 



irrespective of the amplitude of the matter fluctuation 2 . 
This result actually extents to the amplitude of the con- 
vergence three-point function so that one expects to have, 



-0.8 



(15) 



also. In the following the inter-relation between Q and S3 
will be investigated in more details. But it is already obvi- 
ous that a measurement of Q would be useful to constrain 
the matter density parameter of the Universe. 

2.2. The shear three-point function 

The general expression for the shear three-point function 
can be inferred from the Z-space three point function of 
the convergence field. However, because 7 is a 2 compo- 
nent pseudo-vector field there are many possible ways to 
combine shear triplets. The properties of shear patterns 
produced by these triplets may be complex. In order to 
interpret them in a cosmological context it is preferrable 
to focus first on few of them that can be easily inves- 
tigated with analytical calculations. This is the case for 
((7(a7 1 ).7(a3 2 )) l(x')) for a fixed x 2 — Xi separation which 
is a vectorial quantity. This function displays a specific 
pattern which can be viewed as a 2D vector field. 

2.2.1. Computation of ^(^(a:')) 

Let start with a simple configuration, when x\ and x 2 
are pointing in the same location x, for which case we 
have 7 {x 1). 7 (#2) = ^(x). In this case the results make 
sense only when the filtering effects are taken into account. 
Otherwise the result is dominated by the shear behavior 
at arbitrarily small scale and depends on the detailed be- 
havior of the power spectrum and bi-spectrum in the limit 
k — > 00; in a regime that might be totally irrelevant for the 
observations. The shear field is averaged within a smooth- 
ing window, but in order to keep a simple notation, we 
still denote ~f(x) the smoothed shear vector at the win- 
dow location x. For a top- hat window function of radius 
0o, the three-points function can be written as, 

tfixMx')) = 

d 2 i d 2 r w(i do) w(\i + i'\o ) w(i' 0o) p K (i) p K (i') x 

Q K (U') [2cos(20,+r -29 v )e il '-( x '-^u(0 v )+ 



The reduced skewness can however depend on the normal- 
ization of the matter density fluctuations in the intermediate 
regime between the quasi-linear and the nonlinear regime 



cos(20, - 20,0 e i(»+O.(«'-«0 «(0 r+l )l(16) 



where W is the Z-space top-hat window function, 



W(x) 



Ji{x) 



(17) 



where J\(x) is the first Bessel function of the first kind. 

The computation of expression (|l^) is still complicated 
in general, but it simplifies in the large separation limit, 
9o <C \x — x'\. In this case it is possible to express the 
result in terms of the square root of the variance of the 
filtered convergence a K and the convergence-shear corre- 
lation function £ re7 . 

In the large separation limit, in the first term of Eq. 
(|l6|), the values of I' that will most contribute are of the 
order of l/\x — x'\ whereas that of I will be of the order 
of 1/0Q) in the second term the integration in I and V 
factorize out (in this configuration W(\V + l\) does not 
play any role) and both I and V are of the order of l/0o- 
As a consequence the first term of ( |l6| ) is of the order of 
°k £k 7 whereas the second is of the order of £ 2 and is 
thus smaller^ Furthermore, as in the first term V <C I, it 
is possible to expand all quantities in I' /I, 



cos(20 J+r - 20 r ) « 1 - ^ sm 2 (0j - V ) 



(18) 



W(\l + l'\6 ) « W(10 O ) + Wo W'iWo)^ + ... (19) 

Once a shape for Q K is given, computing the integral of 
the angle of I is straightforward. In general, following the 



prescription of Scoccimarro & Couchman, Eq. (11), one 
writeQ 



can 



LI' 1. 1' 



l\2 



Q K (l,l') = a(l,l')+b(l,l') ^— + — J +c (M')^r(20) 
from which one gets, 

{^{x^ix')) = 2n J ldld 2 l'P K (l)P K (l')^' ^'-^ x 

c (i,i'Y 



W 2 (l o )W{l'0 o 
W 



a(l,l') 



W(l0 o )W'(W o )W(l'0o)b(l,l') } u{0 v ). 



(21) 



When a(l, V), b(l, I') and c(l, I') are pure numbers as it 
is the case for the quasi-linear regime, the integrals in I 
and I' factorize and lead to an expression of the form, 

< 7 2 (a07(aO> = c 2 i^ X 

d 2 l'P K (l')W(l'0 o ) e"'-^'-") u(0 v ) (22) 



3 this is a very general scheme for such computations. It is 
described in more details in Bernardeau et al. 2001 

4 The coefficients a, b and c that appear in ( |TT[ ) and (^o|) are 
not necessarily identical because of the projection effects. 
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Fig. 1. Shear patterns in the vicinity of one of the pair point either X\ or X2 (positioned at the center of the plot) 
for n = —1.5, n = —1. and n — —0.5 (from left to right) as a function of x' . It illustrates the results of Eq. (|35|). 
Coordinates are arbitrary. 



where \ is a pure number. It can actually be noted that 
C2 1 identifies with the coefficient appearing in the expres- 
sion of the reduce joint cumulant, 



(k 2 (x)k(x')) = C2\a 2 K (k(x)k(x')) 



in the large separation limit (Bernardeau 1996). It has ac- 
tually been already computed in Bernardeau et al. (1997) 
in the quasi-linear regime. It scales like f^ ' 8 and it is 
mainly independent of the spectrum index. 

The last integral in (^2|) is equal to the convergence- 
shear correlation function, 



£ K Jx')= I d 2 l'P K (l')W(l'9 )e il ' x ' u{9 v ) , 



(24) 



and can be computed by noticing that, 
d 2 l'P K (l')W(l'e ) e [l '- {x '- x) u(p v ) = 
d 2 l' 



2d x >d y i 



I' 2 



■P K (l')W(l%)e 



\V . (x' —x) 



(25) 



So, if one defines £^ as the potential two-point correlation 
function, 



I' 2 



P K (l')W(l'6 )e 



\V .X 



dl' 

7 



- 2tt / —P K (!')W(!%)Jo(l'x), 
we get, 

d 2 l'P K (l')W(l'e )e il ' x ' u(6i,) 



(26) 



x f2 _ y a 
2x'y' 



x 



(27) 



Therefore, in the case of a power law spectrum, P(l) ~ V 
one has for x ^ 8o 



^(x) - x" 
so that, 



(28) 
(29) 



which is equal to — (n + 2)/n ^ K (|a;'|). Hence, the 3-point 
function clearly depends on the slope of the power spec- 
trum: 



( 23 ) (7 2 (*h(*')> 



C2 1 (?l U{\d'\ 



cos(26» d ') 
sin(20 d O 



,(30) 



where d! is the distance vector between x' and x. d! = 
x' — x, and 6^ is its angle to the first axis. It is inter- 
esting to notice that, contrary to the convergence field, 
the amplitude of the three-point shear function vanishes 
(in units of the square of the two-point function) when 
n — » — 2. Although this result is obtained in some specific 
limiting configuration it expresses a general trend: when 
n is close to —2 the shear is dominated by very long wave- 
length^, much longer than a;'— a;, so that what is computed 
here is the same as contracted three point functions that 
all vanish for symmetry reasons. 

For n w —1.5 and sources at redshift unity, we know 
from Bernardeau et al. (1997) that in the quasilinear 
regime, 



C2i «36.7/T> 



0.8 



(31) 



a result that can be obtained from Eq. (|T]) with 
a = 10/7, b = 1, c = 4/7. Consequently, observables 
like (j 2 (x)-f(x')}/ (a 2 £, K (\x'\)) would provide alternative 
ways for measuring the cosmic density parameter fi m . 
They do not require mass reconstruction but still require 
some filtering which, for the reasons mentioned in the be- 
ginning we would like to avoid. 

2.2.2. Computation of ((~f(xi).~f(x2)) j(x')) 

In the previous paragraph the calculations were tractable 
without strong hypothesis on the shape of the bispectrum. 

when n — —2 the computation of the variance of k (or 7) 
shows a divergence at I — > when it is computed with the 
Limber approximation. It means that in the limit n — > —2 
the fluctuations of the shear field are dominated by infinitely 
long wave modes, much longer than \x — x'\. Moreover when 
n < — 2 the whole calculation presented here, which is based 
on the small angle approximation, becomes invalid. 
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Here, we explore more generic geometrical cases, so more 
specific assumptions on the bispectrum are necessary to 
carry out analytical computations. We assume it follows 
the prescription usually adopted in the strongly non-linear 
regime, that is the coefficient Q K introduced in Eq. (|l^) is 
constant (but depends on the cosmological model) . We do 
expect this is a valid approximation at the scales we are in- 
terested in (below 5') and in any case the results presented 
in the following do not critically rely on this assumption 
as it is the case for the skewness of the convergence. 

A quantity like ( (7(0:1). 7(332)) l{x')) is expected to 
behave as (7 2 (33) 7(33')) when \x\ — x 2 \ <C \x' — {x\ + 
x 2 )/2\. In this limit its expression is given by, 

((7(xi). 7 (x 2 )) 7 (x')> = J d 2 ld 2 l'P4l)P K (l')x 

cos(26 l+v - 26 v )u{6 v ) + { Xl «-» x 2 } (32) 

Here the filtering effects can be ignored since all points 
are taken at finite distance. 

In case of a simple bispectrum - with only a non-zero 
monopole term and b = c = - and with the help of the 
expansion (|l8| ) the integrals over the angle 61 and 0j/ in 
(|32j) can be computed explicitly. It leads to an expression 
of the form, 



((■j( Xl ).j(x 2 )) j(x')) = C21 x 

6,(1*2 -Sll) 6c 7 (|d'|) 



cos(20 d O 
sin(20 d /) 



(33) 



where d! = x 1 — (x\ + x 2 )/2. It generalizes the result of 
Eq. ||). 

This result can be simply interpreted: an excess of 
shear at a given position is most likely associated with 
a mass overdensity, so that the shear at finite distance 
from the pair points is preferentially tangential. This non- 
zero three point function is therefore simply associated 
with the usual skewness. It directly comes from the rela- 
tive excess of overdensity regions compared to low density 
areas. 

The situation becomes more complex when the as- 
sumption \xi — x 2 \ <C \x' — (xi + x 2 )/2\ is dropped. 
Another restricting case can however be investigated, 
when \x' — X\\ <C \x\ — x 2 \. In this case the two dom- 
inant contributions are 



£ R (\x' - xx\ 



+ 



(Ctt (I*' 



Xi 



) -Crr(|aj' -SCil)) 



cos(46> x '_ Xl ) 
sm(A0 x >- mi ) 



(35) 



if the x 2 — Xi vector is along the first direction. In this 
expression we have decomposed the shear two-point corre- 
lation function into the tangential part £ tt (i- e - the corre- 
lation function of the shear components along the x 2 — Xi 
direction) and the radial part £ rr . The relative importance 
of the 2 terms depend on the power spectrum shape, e.g. 



Q-2 



£tt-6r (n + 2)(n + 4) 



6, 



(n — 2)n 



(36) 



This shape is actually a direct transcription of the behav- 
ior of the shear two-point functions (for a component 71 
of 7(031) along a fixed direction), 



(71(^1) = 



£ K (\x' - xi| 



(t[ tt (\x' - Xl \) - Z TI (\x' - Xl \)) 



cos(46 x ,- Xl ) 
sin (4d x x J 



.(37) 



((7(xi).7(x 2 )) 7(33')) = / d 2 ld 2 l'P K (l)P K (l')Q K (l,l') x 
s u.(«a-«i)+u'.(«'-«0 C08 (2Q l+v _ 29i) u{6 v )+ 
e iI -(»'-»0+«'.(»'-««a) COB (2fl | _2ft/)«(e|/ +I )l (34) 

Then similar calculations can be performed, again assum- 
ing the bispectrum simply factorizes in terms of the power 
spectrum as in the nonlinear regime. It leads to the ex- 
pressions, 



((7(xi).7(a; 2 )) -y(x')) = e 2 i £« 7 (| 



332-331 



Nonetheless these patterns are indicative of the physi- 
cal effects one might want to look for. The pseudo- vector 
field structure of the three-point function displays a pat- 
tern which corresponds to the superposition of a uniform 
field and a specific mass-dipole contribution. The relative 
weight of the two terms depends on the power law index. 
When n is close to -2, structures are mostly dominated by 
long- wavelength modes and the shear patterns are aligned 
(in other words £tt and £ 1T are equal). In contrasts when 
n is close to 0, structures are given by point-mass Poisson 
distribution (and £tt and £ rr are opposite to each-other). 

Plots on Fig. ^show the shear patterns near pair points 
for different power law indices. One can see that for n < 
— 1 the shear pattern around the points is mostly uniform. 

2.3. Semi-analytical results 

The exact analytical results presented in the previous 
paragraphs only correspond to simple cases. They pro- 
vide insightful descriptions of the behavior of the shear 
three-points function but have a limited practical interest 
if their exact validity domain is unknown. More general 
results can be obtained, but only for specific cosmological 
models through numerical computation and ray-tracing 
simulations. 

Fig. ^ shows the shear pattern expected for a power 
law spectrum of index n = — 1 and for a matter three- 
point function shape given by the non-linear regime. The 
patterns observed around the pair points (located at po- 
sitions [-1,0] and [1,0]) are rather complicated. In partic- 
ular the circular shape expected around the pair points 
(according to the results of sect. 2.2.2) are observed only 
at very large distances (more that 5 times the pair sep- 
aration). At separation comparable to the pair distance 
the patterns increase complexity, with a substantial part 
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Fig. 2. The shear pattern obtained by a numerical inte- 
gration in case of a power spectrum index n = — 1 and for 
a three-point function that obeys the hierarchical ansatz, 
Q K = constant in Eq. (fl2|). The pair points are at position 
(—1, 0) and (1, 0) in both panels. The only difference in the 
two panels is the scale at which the pattern is drawn. 



of the area with no significant correlation. The most strik- 
ing feature is the quasi-uniform shear orientation along the 
segment joining the two-points. This effect which is indeed 
expressed in Eq. (|3^), where a uniform component is ex- 
plicitly predicted, clearly strengthens in between the pair 
points. Such a feature might appear somewhat surprising 
but a close inspection of synthetic shear maps (Fig. ||) 
indeed reveals a lot of highly contrasted clumpy regions 
surrounded by strong coherent shear patterns primarily 
oriented transversely to directions between lumps. 

The central pattern is the strongest and the most typ- 
ical feature of the three-point correlation map and should 
be the easiest detectable one in the data. The previous an- 




Fig. 3. Example of shear map obtained from a numerical 
simulation (open CDM model, for an 5x5 square degree 
survey, van Waerbeke et al. 1999) 



alytical results suggest that such a structure is expected to 
hold for power law indices between —2 and —1 and should 
then be robust enough to be used as a detection tool of 
non-gaussianities. More detailed analysis performed along 
this idea are presented in the following. 



3. Comparison with numerical simulations 

The numerical simulations we use are described in Jain 
et al. (2000). The cosmological model is an open Universe 
{fl m — 0.3, £1a = 0) with a Cold Dark Matter power 
spectrum (r = 0.21) and a normalization erg = 0.85. The 
sources are located at redshift unity and the simulation 
area covers about 11 square degrees with a resolution of 
0.1 arcmin. 

The shear patterns for pair points at increasing sepa- 
ration is shown on (Figure ^) . A visual inspection of their 
morphology and strength confirms that the uniform shear 
pattern within an ellipse that encompasses the pair points 
(as described in the next section) is likely an optimum 
way to extract non-gaussian signal. When the separation 
is small, the overall circular shear pattern is clearly visible. 
When the separation increases the shear appears uniform 
in the neighborhood of the segment joining the pair points, 
and is mostly radial at finite distance. We have already 
seen that these results might somewhat be dependent on 
the power spectrum index. For this simulation the index 
varies from —1.3 to —1 and it is thus natural that the 
patterns look like those obtained in case of a power-law 
model n = — 1. 
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Fig. 4. The shear pattern of ((■y(xi).'y(x2))'~f(x')) mea- 
sured in ray-tracing simulations (Jain et al., 2000) for in- 
creasing pair separation \x2 — x±\. The separations are, 
from top to bottom and left to right, 2, 4, 6, and 8 in plot 
units (1 unit corresponds to about 10"). The pair points 
are along the horizontal axis. 

4. Improved measurement strategies 

In this section we compare the measurements made in 
mock catalogs that mimic a large number of observational 
effects with different input models. We use these results 
to develop different survey strategies adapted to real data 
set. 

4.1. Mock catalogues 

The mock catalogues are generated from simulated sky 
images following the procedure described in Erben et al. 
(2001). The only difference here is that the galaxies are 
lensed according to a realistic cosmic shear signal, using 
ray-tracing simulations (Jain et al. 2000) instead of having 
a constant shear amplitude as in Erben et al. 2001. The 
galaxies are analyzed exactly in the same way as real data, 
following the procedure described in van Waerbeke et al. 
(2000, 2001b). In particular, the mock catalogues contain 
the main three features encountered in the actual surveys: 

— galaxy intrinsic shape fluctuations; 

— masks; 

— noise from galaxy shape measurements and systemat- 
ics from PSF corrections , 

and the simulated images reproduce observational con- 
ditions matching our data set (PSF anisotropy, limiting 
magnitude, luminosity functions, galaxy and star num- 
ber densities, intrinsic ellipticity...). The PSF anisotropy 
is even larger than the one in our data set, as shown in 




<e5,„> ,/3 



Fig. 5. Histogram of the Point Spread Function 
anisotropy of the stars in the simulated images (solid line) 
and in the VIRMOS-DESCART survey (dashed line). 

Figure [|, but it is uniform [j. We used two ray-tracing 
simulations from Jain et al. (2000): one is OCDM, as 
described in Section 3, and the other is a rCDM with 
r = 0.21 and O m = 1. For each simulation we produced 11 
square degrees of simulated sky images containing roughly 
30 galaxies per arcmin 2 , with a pixel size of 0.2 arcsec. The 
galaxies are lensed following the ray-tracing shear before 
generating the image |^. In addition to the galaxy mock 
catalogues, we produced a reference catalogue containing 
only the cosmic shear values. We can therefore study sep- 
arately the effect of masks, the ellipticity Poisson noise, 
real noise and systematics, when compared to the refer- 
ence catalogue results. 

4.2. Detection strategy 

The adopted strategy is to measure the two component 
quantity 

U\*i-*,\)= I /^<(7(zi).7(*2))7(z')} (38) 

./ellipse ^ellipse 

for different separations. The bar means we perform the 
average of ((■j(x 1 ).'y(x2)) ~f(x')) within the area defined 
by Veiiipse- The ellipse properties are defined to cover an 
area that encompasses the pair points in which the ori- 
entation of the shear pattern is expected to be uniform. 
To be more specific, after some trials we have obtained 
a good balance between the signal intensity and the shot 
noise amplitude when the ellipse has a 0.4 eccentricity, the 
pair points being its foci (see on Fig. ^). 

The computation time can be reduced by building first 
a Delaunay triangulation (Delaunay, 1934; see also van de 
Weygaert, 1991) of the survey so that the neighbor list 
of any galaxy in the survey can be easily (and quickly) 

6 However, we do NOT assume a uniform PSF for the PSF 
correction, as for the real data, we fit the PSF with a 2D second 
order polynomial 

7 That is the image of the galaxies are not lensed. We first 
lens the catalogue of source galaxies and then use it to generate 
the sky images 
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Fig. 6. Results for the measured value of the two-point function and the reduced three-point function (as defined in 
Eq. [58]) as a function of the pair distance for an Open CDM model (thick dashed lines) a ACDM (thick dot-dashed 
lines) and a rCDM model (thick solid lines). Results obtained in an Open CDM models with masks is shown as thick 
dotted lines. Galaxies in the catalogs have no intrinsic ellipticities. The measurements are made in bins of width 10" . 
The thin lines of the right panel corresponds to 200 (n + 2)/(— n), 150 (n + 2)/ (— n) and 60 (n + 2)/ (— n) for the open, 
A and r CDM case respectively, where n is the power spectrum index measured from the slope of the corresponding 
two-point correlation functions. 



constructed. The quantity £3(1x1 — a^l) is a pseudo- vector, 
but its second component, e.g. the one corresponding to 
elongation at 45 degrees to x\ — x 2 , vanishes in average 
for symmetry reasons. Therefore we only have to compute 
one component, which we define as, 



£3 ( d ij) 



y.,,1.- "V "'./ "> (ei- e j) e 

Eijfc W * W 3 w k 



(39) 



where e, are the galaxy ellipticities, ejj?^ is the tangential 
(e.g. first) component, with respect to the segment Xj—Xi, 
of the ellipticity of the galaxy of index k, to, are weights 
associated with each galaxy. The sums are made under the 
following constraints: 

— dik > d m i n , djk > d m i n so that close pairs are ex- 
cluded to avoid spurious small angular scale signals 
(van Waerbeke et al. 2000) and 

— dik + djk < 1.1 dij so that galaxy fc is within the afore 
defined ellipse and the sum is splitted in bins according 
to the d^ distance. 

We then plot the reduced three point function, that is 
£ 3 in units of [6(<2y)] 2 , 



6 ((k 



^EijWjWj (ej.ej) 



(40) 



The noise is estimated by computing the r.m.s. of the 
estimator in small individual bins (for which the noise 
dominates the signal), and then rebinned in the final, 
larger, bins (as described in Pen et al., 2002). The noise 
of the reduced three-point function is computed assuming 
that the noise of the second and the third moments are un- 
correlated. This might not be exact but the consequence 



of this assumption is insignificant since the main contribu- 
tion to the noise of the ratio is in general the three-point 
moment. 



4.3. Behavior of the reduced three-point function 

In Fig. ^ we show the behavior of the reduced three-point 
function in the simulations. These simulations do not con- 
tain any noise, except that the shear field is measured in 
discrete randomly placed points. It is found to be a slowly 
growing function of scale. This might first appear surpris- 
ing (this ratio is expected to be constant for a constant Q) 
but this result should be examined in the light of equa- 
tion ( |33| ) and the subsequent comments: at small scale n 
decreases, therefore it is expected that the reduced three- 
point function also decreases. Indeed, as demonstrated on 
the right panel of Fig. |[ the reduced ratio behaves approx- 
imately like (n + 2)/n at least for scale above 2' with a 
multiplicative coefficient (200, 150 and 60 for respectively 
an open, A and r CDM model) that roughly correspond 
to the values of the convergence skewness S3 at small scale 
in those models (Hui 1999). 

On this plot the effects of masks on these quantities are 
shown as dotted lines. As expected its effect on the two- 
point function is null. Its effect on the three-point function 
as measured here is not totally absent (different configu- 
rations taken into account in Eq.(|3^) may have different 
weights when masks are taking into account). It appears 
however that it has a negligible effect compared to the 
other sources of noise. 

Fig. [?] also shows the amplitude of the cosmic variance 
to be expected in an 11 square degree survey for the re- 
duced three-point function and in case of an Open CDM 
model (dotted dashed lines). These error bars have been 
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obtained from the results obtained in 7 realizations of the 
same model. They show that 30% fluctuation is to be ex- 
pected in the signal. It is to be noted also that the error 
bars are correlated in the different bins. 

4.4. Optimization and noise effects 

A proper choice of weights Wi is essential for getting a good 
signal to noise ratio. We have found that the introduction 
of a cut-off in the galaxy ellipticity distribution, 



exp 



2 el 



(41) 



improves the signal detection but may also affect its am- 
plitude. In figure ^| a cut-off is introduced with eo = 0.5 
in order to improve the S/N ratio. If the measured shear 
is simply the sum of the intrinsic shear and some noise, 



(42) 



then the cut in e translates in a cut in 7 that can be simply 
described. The result depends on the shape of the intrinsic 
ellipticity PDF, P(e), 



<v cut = 
/ 1 



/d 2 e ( 7(a;4 ) + ei )F( e ) exp(-(^) 
Jd^eP(e) exp(-^ 



(43) 



If the cut-off value eo is large enough compared to the 
typical excursion of cosmic shear values, then this equation 
can be linearized in 7 and one finds that, 



7r* = /c7«>+ 



with 



fc = 



fede(l- 



2rr 



P(e) exp 



2ei 



JedeP(e) exp (-^) 



(44) 



(45) 



For the adopted ellipticity distribution and for eo = 
0.5, the value of f c is about 0.8. Moreover the cut-off is 
large enough to have no impact on the non-gaussian prop- 
erties of the shear field (e.g. sub leading order terms in 
Eq.(f44|) have a negligible effect). 

These numerical investigations are useful for exploring 
the effects of various noises. As expected the intrinsic el- 
lipticities increase the noise, but do not bias the result. 
For 11 square-degrees surveys the signal to noise for the 
detection in a rCDM model is still comfortable (see Fig. 
I?]). Masking effects are also rather mild. They bias very 
weakly the results (see Fig. |^) and only slightly deterio- 
rate the signal to noise ratio when intrinsic ellipticities are 
taken into account (because the number of available pairs 
or triplets is smaller). 

Moreover, when similar measurements are made in 
mock catalogues that incorporate realistic noises we de- 
scribed above, we recover the expected signal with still a 
reasonable signal to noise ratio: it is still larger than 3 in 3 
independent bins (see Fig. ^). Note that the S/N ratio for 



1000 r 



100 r 




d 12 / arcmin 

Fig. 7. The reduced three-point function as in previous 
figure for different noise types in the catalogs. Error-bars 
for the open CDM model (dotted-dashed line) correspond 
to a cosmic variance estimation. The plot also compares 
the case of exact ellipticities for a rCDM model (solid 
lines) with the case where galaxies have intrinsic elliptici- 
ties (dotted lines, the measurement is made in 4 indepen- 
dent bins) and the case where noises in the shape determi- 
nation are included (after PSF corrections, etc., dashed 
lines, 4 bins). Error bars correspond here to measurement 
errors, not to the cosmic variance, and are therefore un- 
correlated. 

the reduced skewness is not only sensitive to the ampli- 
tude of the three point function, but also to the amplitude 
of the 2-point function. On the lower curves of Fig. [?], only 
the noise coming from the three-point function measure- 
ment is shown. Its amplitudeis roughly independent on 
the cosmological model (it scales basically like the inverse 
square root of the number of triplets available). 

5. Conclusion 

We have presented various geometrical patterns that the 
shear three-point correlation function is expected to ex- 
hibit. Its dependence with the cosmological parameters is 
found to be similar to that of the skewness of the con- 
vergence field, opening an alternative way to break the 
degeneracy between the amplitude of the density fluctua- 
tions, eg, and the density parameter of the Universe f2 m . 
However, the reduced three-point function of the shear 
is more dependent on the power spectrum index than 
the skewness. In particular it is expected to vanish when 
the index gets close to —2. Numerical investigations have 
nonetheless proved the shear patterns to be robust enough 
to provide a solid ground for the detection of non-gaussian 
properties in cosmic shear fields. 

We proposed a detection strategy that has been tested 
in mock catalogues that include realistic noise structures 
such as residual systematics and PSF anisotropy as seen 
in the real data. The quality of the PSF correction is al- 
ways good enough to provide an accurate measurement 
of the shear three-point function. Since the mock cata- 
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logues were designed to reproduce the characteristics of 
the current VIRMOS-DESCART lensing survey, we con- 
clude that non-gaussian signal should be detectable in this 
data set. The result of our investigations is presented in 
another paper (Bernardeau et al. 2002). 

Beyond the detection, the scientific exploitation of the 
3-point function for cosmology also depends on our ability 
to overcome other important difficulties. For instance, we 
found the cosmic variance amplitude to be of the order of 
30% of the signal on the reduced three-point function, in 
agreement with previous studies made for the convergence. 
Other issues regarding source clustering, source redshift 
uncertainties, and intrinsic alignment of galaxies have not 
been considered here. 

The measurement of non-gaussian signatures in lens- 
ing surveys is of course of great interest because it pro- 
vides an independent measure of the mean mass density 
of the Universe in addition to test the gravitational in- 
stability paradigm which lead to large scale structures. 
It is likely that the analysis of cosmological non-gaussian 
signatures will be one of the major and most promising 
goals of emerging dedicated lensing surveys^] that take ad- 
vantage of panoramic CCD cameras of the MEGACAM 
generation (Boulade et al 2000). 
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